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We consider the production of charm by real and virtual photons. Special attention 
1 is paid to diffractive charm production, which provides information on the gluonic content of 

the Pomeron. Our calculations are based on the gluon distributions of the CKMT-model, 
which is shown to lead to agreement with the data on open charm production in deep 
inelastic scattering. We compare predictions for diffractive charm production of different 
^ \ models for the distribution of gluons in the Pomeron. Experiments at HERA should be 

able to discriminate between them. Predictions for beauty production in diffractive and 
non-diffractive interactions of photons are also given. 

o 

On 

Cii! I. INTRODUCTION 

The production of heavy quarks in deep inelastic electron scattering from a proton has received 
increasing attention recently, since this reaction is seen as a tool to probe the gluon distribution in the 
nucleon. Due to the large mass of the quark, the reaction is believed to be driven by a perturbative 
mechanism, photon-gluon fusion, and is therefore sensitive to the nucleonic gluon density. We consider here 
in a consistent way the production of charm in deep inelastic scattering (DIS) as well as in hard diffractive 
scattering. The latter process, while driven by the same basic reaction mechanism, probes a different gluon 
density. As the diffractive processes are mostly seen as mediated by the Pomeron exchange, the relevant 
quantity in diffractive production is then the gluon distribution in the Pomeron; a review can be found in 
Ref. @ . 

The distribution of gluons in the proton is known now comparatively well from deep inclatic scattering 
for x > 10~ 2 , but for smaller values ofx information on g p (x, Q 2 ) is very limited. A comparison of theoretical 
predictions on charm production in the HERA energy range with experimental data allows one to test small- 
x behaviour of gluonic distributions. 

Diffractive production of hadrons in deep inelastic scattering has been observed in experiments at 
HERA [|]||. The gluonic content of the Pomeron is poorly known. Studies of the Q 2 -dependence of the 
Pomeron structure function lead to the conclusion that the distribution of gluons in the Pomeron is 
hard and that they carry the main part of the Pomeron momentum. However there are big differences 
between gluon distribution functions gp{x,Q 2 ) in different models. 

It is the purpose of this paper to present a consistent approach to open charm photo- and electro- 
production as well as diffractive charm production. The gluon distributions we use are based on the model 
of Ref. HJ, hereafter referred to as CKMT. No new parameters are introduced into our calculations. In this 
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approach there is only one basic Pomeron, which contributes in different ways to both "soft" and "hard" 
mechanisms. This can be represented in terms of an effective or "dressed" Pomeron with a Q 2 -dependent 
intercept. It was the basis of a very compact parametrization by CKMT for the structure function F2 and 
distributions of quarks and gluons in the proton, which was used as an initial condition for QCD evolution 
||. Also diffractive vector meson production (p, (f>, ip) by real and virtual photons [|o| could be described 
very well by this model without additional parameters. 

There have been a variety of investigations on aspects of charm production in deep inelastic and 
diffractive scattering, in part going beyond the leading order. In order to achieve a clear and simple 
discussion of these reactions, we confine ourselves to the leading order terms in the cross sections. This 
more phenomenological approach seems reasonable in view of the large sensitivity of perturbative QCD 
calculations to the heavy quark mass and the present uncertainty in the data. 

In Section 2 we shall first consider inclusive charm photo- and electroproduction on the proton. We 
shall formulate a method of calculation of the cross sections of heavy quark production for arbitrary values 
of Q 2 . It will be shown that the gluon distribution proposed by CKMT gives a very good description of 
the existing data on charm production. We therefore proceed in Section 3 to calculate charm production 
in diffractive scattering in the context of the CKMT model and compare to other models. In Section 4 we 
present predictions for inelastic and diffractive beauty production by real and virtual photons. 



II. CHARM PHOTO- AND ELECTROPRODUCTION 



In this section we consider the open charm production by real or virtual photons on a proton. In 
the past several prescriptions have been proposed to incorporate the contribution of charm into F2(x,Q 2 ). 
A widely used method |ll]-|l4j was to generate the charm quarks dynamically, starting with no intrinsic 



charm below some threshold Q 2 h ~ m 2 and to produce charm quarks (considered as massless) through QCD 
evolution. However this results into a too large charm contribution in the threshold region. We therefore 
choose a different approach. For moderate Q 2 , the fact that the mass of the charm quark is comparatively 
large, m c ~ 1.5 GeV 3> Aqcd, makes it possible to apply perturbation theory. It was shown in Ref. [l5| 
that therefore more realistic predictions for charm production can be obtained from the photon-gluon fusion 
diagram of Fig. la. This is true in a broad region of Q 2 . Next to leading order (NLO) calculations in QCD 
perturbation theory have been carried out [p"5|-p2"| for heavy quark production. These higher order terms 
lead to no new qualitative feature and can be incorporated by adjusting the parameters of the lowest order 
calculation Jlq ] . In view of this perturbative stability and since there are uncertainties in the charm quark 
mass and the factorization scale, we prefer to simply work in the leading order (LO), where the calculations 
are very transparent. We shall show that by properly chosing the charm quark mass and factorization 
scale one obtains a good description of the data for charm production in DIS. For Q 2 3> m 2 , on the other 
hand, higher order diagrams of the QCD perturbation theory should be resummed [2j|. Since mass effects 
can then be neglected, the charm quarks will be produced dynamically through QCD evolution, where the 
charmed quark is considered as massless. We obtain the charm quark distribution needed as input for this 
evolution through the photon-gluon fusion mechanism. 

At moderate Q 2 , the main contribution for large centre-of-mass energies a/s is expected to come 
from the gluon fusion process, shown in Fig. la, where a gluon from the proton interacts with the photon 
and produces a charm anti-charm quark pair. We will discuss this mechanism first. Contributions from 
"resolved" charm production, Fig. lb, are less important and will be dealt with later. 

In QCD perturbation theory the cross section for the process ^*'p — > ccX can be written as a 
convolution of the gluon distribution g p (z,p 2 ) and the partonic cross section of the photon-gluon fusion 
process (7(7(7 — > cc), 

f 1 Q 2 

a{-f*P -> ccX) = / dza 7g ^ c5 (x,z,—,ml)g p (z,p 2 ) . (1) 
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The threshold condition for production af a charm-anticharm pair leads to a lower bound for the fraction 
z of the proton momentum carried by the gluon, z m i n = ax, where a = 1 + Am 2 /Q 2 and x is the Bjorken 
variable. For the partonic cross section of charm production by photon-gluon fusion, one finds [^3 24 with 
the Hand convention for the flux 

a (x z m 2 ) - ^ aEM 2xc 2 as{ ^ F) C ( - (2) 

where /x/ is the factorization scale and e c the charge of the charm quark in units of e. In leading order the 
coefficient function C is 

C(C, r) = \ [C 2 + (1 - C) 2 + 4C(1 - 3C)r - 8C V] In \±1 
2 1 — v 

+|[-1 + 8C(1-C)-4C(1-C)r] , (3) 

with 

^ 2 = 1-^. (4) 

In the limit Q 2 — Eq. ([!]) yields the photoproduction cross section. The contribution of charm production 
to the proton structure function is given by 

Ff(x, Q 2 ) = Q ^ ~ X) <t( 7 *p -> ccX) . (5) 



For the gluon distribution, a crucial ingredient in the calculation of the cross section, we use the 
CKMT-parametrization of Ref . || . It gives a good description of the HERA data for the proton structure 
function F2(x,Q 2 ) at all Q 2 and for diffractive vector meson production p,[[o|]. We briefly summarize here 
only the physics background; details can be found in the Appendix and in Ref. H. 



Experimental studies of small x DIS at HERA P5|]26[ had shown a fast increase of the proton structure 
function Fq.{x, Q 2 ) as x — > 0, which was considered by some authors as an evidence for a "hard" Pomeron. 
This Pomeron has an intercept Op(0) substantially above unity, contrary to the "soft" Pomeron, which is 
observed in high-energy hadronic interactions and photoabsorption, which has an intercept only slightly 
above 1: ap(0) = A + 1 w 1.08. It was argued however by CKMT in Ref. B that this ap(0) extracted 
from high-energy behaviour of hadronic total cross sections is not the actual intercept of the Pomeron itself, 
but an effective value that incorporates large effects of Pomeron rescattering (multi-Pomeron cuts) in soft 
processes. The actual "bare" value of the pole intercept extracted from the analysis of many features of 
hadronic interactions, taking into account multi-Pomeron processes, was found to be substantially higher, 
corresponding to A w 0.2 (27|. At large Q 2 in DIS the contributions from rescatterings are much smaller than 
in hadronic interactions (or ^yp) and the "bare" intercept determines the behaviour of structure functions 
at Q 2 > 1 GeV 2 . In this approach there is thus only one Pomeron, which contributes to both "soft" and 
"hard" processes. This effective or "dressed" Pomeron then has a Q 2 -dependent intercept. It was used by 
CKMT to provide a compact parametrization of the structure function F2 and the distributions of quarks 
and gluons in the nucleon. These distributions are used as initial condition for the QCD evolution 

The gluon distribution resulting from this model is 

xg p {x,Q 2 ) = C g x-^\l-xY^ 2 ^ , (6) 

with the effective Pomeron intercept determined by 
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a «5 2 »- A <°»( 1 + i?T*) ' < 7 > 

For low x the gluon density in this model exhibits the x dependence which is characteristic for the Pomeron 
exchange. The behaviour at x ~ 1 has been obtained from counting rules for the sea quarks which are 
a factor 1 — x softer than the gluons. Finally, the factor C g follows from the momentum sum rule. This 
parametrization is valid up to Q 2 = 5 GeV 2 . We actually use it only up to 2 GeV 2 and for higher values 
QCD evolution is applied where the parametrization provides the initial condition. We take the running 
coupling constant with A = 0.20 GeV and four flavours. 

The results for charm production are very sensitive to small changes in the value of the charm mass. 
We will use charm photoproduction to fix the value of this parameter. In Eqs.(|l|) and (0) one encounters 
the factorization scale /i/. The NLO results in Ref. jl5| showed that the scale /z 2 = 4m^ yields the best 
stability of the pcrturbative calculations and we will use this now for photoproduction; other choices will 
be discussed below in connection with electroproduction. Since this factorization scale is above the value 
of Q 2 for which the parametrization is valid, we use the leading order QCD evolution equations for three 
flavours to obtain the gluon distribution at fi 2 . 

We calculated open charm production corresponding to the diagram of Fig. la with the gluon distri- 
bution function as given in (^|). The results for the total inclusive cross section are shown for four different 
values of the charm mass in Fig. 2. One can see that relatively small changes in the charm mass lead to 
large differences in the cross section. For a charm quark mass m c — 1.4 GeV our curves are close to the 
experimental data and look similar to the results of the NLO calculation jnj^2) ■ I n the following we shall 
use this value of the charm quark mass. 

For very large energies y/s the resolved contribution, shown in Fig. lb, starts to be important. In this 
production mechanism a gluon from the photon interacts with a gluon from the proton to produce a charm 
anti-charm pair. To obtain the contribution from the resolved production, we therefore need the gluon 
distribution in the photon, g^{x,Q 2 ). This distribution is poorly known at present. For an estimate of the 
resolved contribution, we use a distribution obtained by the same method as for the determination of the 
gluon distribution of the proton in Ref. |9| and thus introduce no new parameters. At small x factorization 
takes place in the Regge pole model. Therefore the gluon distribution in the photon at low a; is proportional 
to that of the proton; for this proportionality factor we find 

e 7 = Um frfr'Q*) = lim 41 « 0.003 . (8) 



This factorization is approximately valid also when Pomeron cuts are taken into account. The behaviour in 
the limit x — > 1 is different for the photon and proton and can be found from counting rules. So we write 
g-y(x, Q 2 ) in the form 

x 9l (x,Q 2 ) = e ffi^ff = eJC g x^ 2 Hl - x)^^ . (9) 



The contribution of the diagram in Fig. lb is shown in Fig. 2. It is negligible at energies up to 
yfs ~ 10 GeV, but increases with energy faster than the dominant photon-gluon fusion contribution of 
Fig. la. This contribution is of the order of 10% already at HERA energies. For other models of <? 7 (x, Q 2 ) 
this contribution is larger p2|. It is also sensitive to the value of the factorization scale and increases for 
smaller values of this scale. 

We now turn to electroproduction. First we will consider the factorization scale 4m 2 , as in photopro- 
duction. The production of charm quarks give a large contribution to the DIS structure function F 2 (x, Q 2 ) 
at small x and large Q 2 m 2 , where they can be considered in the framework of QCD evolution equation 
as massless. However at Q 2 ~ 10 GeV 2 or less the charm mass is very important and the charm quarks 
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cannot be treated in the same way as the light quarks |15| . We therefore deal with electroproduction in 
two different ways, depending on the value of Q 2 . For values of Q 2 less than a certain value Q 2 , we use 
the perturbative approach for the gluon fusion process with massive charm quarks as given by Eq.([lj). To 
obtain the gluon distribution at the factorization scale Am 2 , we proceed as in the photoproduction case; 
we start the QCD evolution for three massless flavours from Q\ — 2 GeV 2 . With this distribution, the 
charm contribution is directly obtained from Eq.(||). For virtualities larger than Q 2 the charm quarks are 
produced through massless QCD evolution with four flavours. For the input distributions, we use the light 
quark distributions as obtained above from three flavour massless QCD evolution up to Q 2 . The charm 
input distribution is generated by photon-gluon fusion at Q 2 . 

We determine the transition value Q 2 from the perturbative approach with massive charmed quarks 
to the massless QCD evolution by demanding that this procedure creates an F% c below and above Q 2 
with a smooth derivative with respect to Q 2 at this point. We found that the values of Q 2 yielding a 
smooth transition in the region of small x vary from 30 GeV 2 to 100 GeV 2 . In Fig. 3a we show the charm 
contribution to the proton structure function, F 2 (x,Q 2 ), as given by Eq.(||) as a function of the virtuality 
Q 2 at different values of x; the transition from massive to massless treatment of the charm quarks is made 
at 50 GeV 2 . We also show the logarithmic derivatives in Fig. 3b to see their discontinuities at Q 2 . The 
change of A in the running coupling constant at the charm threshold has been taken into account. It follows 
from Figs. 3a and 3b that for very small x (x < 10 -3 ) the transition from one regime to the other is very 
smooth and for Q 2 > 50 GeV 2 charmed quarks can safely be considered as massless in the QCD evolution 
equations. This is in agreement with results obtained in Ref. Jl9] and also with Ref. fl5fl , where it has been 
stated that for W 2 — Q 2 (l — x)/x < 10 6 GeV 2 the gluon fusion model should be applied. For x ~ 0.1 the 
mass effects are important up to much larger values of Q 2 . Thus, from Fig. 3a we can conclude that for the 
small x and Q 2 < 500 GeV 2 the difference between the prediction of the gluon-photon fusion diagram and 
its QCD-evolved contribution is small. 

So far, we have used the factorization scale 4m 2 . Another natural candidate for electroproduction is 
4m 2 + Q 2 . To examine this possibility, we show in Fig. 3c and 3d the analogous results for this scale. The 
features are quite different. Now for large x the discontinuity is small, but for small x it is large, just the 
reverse of what one finds for the constant factorization scale. This suggests that one could use a different 
factorization scale and/or different transition value Q 2 for different kinematical regimes to obtain a smooth 
transition from massive to massless quarks. 

Finally, in Fig. 4 we make a comparison of our predictions for charm photo- and electroproduction 
with the data. We show results as a function of W — y/s for the factorization scales 4m 2 and 4m 2 + Q 2 ; 
for the latter we only show predictions at high Q 2 where the difference to the fixed scale is large. Note 
that the cross sections for the different values of Q 2 in Fig. 4 are rescaled. The photoproduction data at 
high energies are from ZEUS and HI [ p8|j29|: the other experiments are listed in Ref. |2g) l. The low energy 
data for electroproduction are from Ref. |30| , |3~l| and the high energy data from HI [32[ . The latter data 
correspond to Q 2 values slightly different from those of the low energy data and of our calculations. We see 
that all data are rather well described. The high Q 2 data indicate a preference for the fixed factorization 
scale of 4m 2 : the charm production at different (but not very large) Q 2 depends on the gluon distribution 
at a fixed scale. Therefore we will use this factorization scale also for diffractive charm production. 

The deviation of the HI results at Q 2 — 12 GeV 2 from our theoretical prediction is partly explained 
by the fact that for our prediction we took the same Q 2 as for the low energy data, i.e. Q 2 — 13.9 GeV 2 . 
It means that our predictions are here slightly too small. Nevertheless, the discrepancy remains significant 
and is difficult to explain in view of the fact that for other values of Q 2 the agreement between theory and 
experiment is reasonable and that the cross sections have only a weak Q 2 dependence in this region. Except 
for these HI points, our comparison confirms the photon-gluon fusion mechanism in combination with the 
CKMT gluon distribution at the factorization scale 4m 2 . 

In comparing our results to the data, it should be noted that the "resolved" diagrams of the type 
shown in Fig. lb can contribute also for highly virtual photons. Our estimates indicate that this contribution 
at present energies is rather small — it is less than 10% of the main diagram of Fig. la. However in the same 
way as for real photons the contribution of Fig. lb increases with energy faster than for Fig. la. Thus it will 
dominate at superhigh energies where it corresponds to charm production in the central rapidity region. 
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III. HARD DIFFRACTIVE SCATTERING AND CHARM PRODUCTION 



In deep inelastic scattering, there are events where particles are produced only in the fragmentation 
regions of the initial particles. We consider here the case that the proton remains intact, the so called "single 
diffractive scattering" . This corresponds to a large "rapidity gap" between the diffractively produced states 
(with invariant mass Mx) and the recoil proton. This observation is readily explained in the Reggeon theory 
through the exchange of the Pomeron (Fig. 5a). 

In this section, we discuss the diffractive cross section in the Pomeron exchange model p^ ] and 
specifically consider diffractive charm production. As in the total inclusive charm production discussed 
above, gluon fusion is again the dominant production mechanism for charm. Therefore, diffractive charm 
production is sensitive to the gluon distribution in the Pomeron, which is believed to consist mainly of 
gluons. This should be compared to the analogous inclusive process, where mainly the gluon component of 
the proton is probed. We will use the CKMT approach to obtain the parton distributions in the Pomeron 
and compare its results to two other models for the Pomeron parton distributions. 



The differential production cross section can be written as 



,2 



-\F?{x,Q 2 ,x P ,t) , 

,Xp,t)\ J 



1 ~ ^ + oTi I d 7—7y> ITT >F?(x,Q*,x P ,t) , (10) 



dxdQ 2 dx P dt xQ 4 { 2 [1 + Rbif(x, Q 

where we have introduced a diffractive structure function F®(x,Q 2 ,xp,t). Here x,y and Q 2 are the usual 
DIS variables and t is the squared momentum transferred to the proton, t = k 2 = (p — p') 2 . The variable 
xp is defined by 

q-k M 2 X + Q 2 
XP= - W2 T D2 ' 11 

q ■ p W z + <4 Z 

with M\ = (q+k) 2 the squared invariant mass of the diffractively produced particles and W 2 = s = (q+p) 2 , 
the squared CMS energy of the photon-proton system. For the mechanism shown in Fig. 5a, the variable xp 
can be interpreted as the fraction of the proton momentum carried by the Pomeron. The function -Rdif in 
Eq.(flC)|) is the ratio of the longitudinal to transverse part of the cross section. Integrating F®(x,Q 2 ,xp,t) 
over xp and t we obtain the diffractive contribution to the total deep inelastic structure function F^ix, Q 2 ). 

For the exchange of the Pomeron pole, the diffractive structure function F®(x,Q 2 ,xp,t) can be 
factorized into two parts, 

F.?(x,Q 2 ,x P ,t) = f(xp,t)Fp((3,Q 2 ,t) , (12) 

where the variable (3 = Q 2 /(M X + Q 2 ) — x/xp plays the same role as the Bjorken variable x has in DIS: 
the Pomeron momentum fraction carried by the partons in the Pomeron. The first factor represents the 
Pomeron flux from the proton and can be written in the form 

f(x P ,t)^ { -&^ Xp ~ 2 ^ , (13) 

1D7T 

where g^ p denotes the Pomeron-proton coupling. The value of ap (t) in this flux factor should be taken at 



some effective virtuality scale, Q e //i since in the model of Ref. § one has 

a P {t) = l + A(Ql ff )+a'pt , (14) 

where A is given by Eq.(|^). The scale Q 2 ff is a priori not known, but it was argued in Ref. § that it 
should be low, since a hard scale from the top part of Fig. 5b does not get through to the lower part of 



G 



the diagram. From theoretical point of view, values of A(Q 2 ^) = 0.13 to A(Q 2 ^) = 0.24 are possible, 
corresponding to the effective Pomeron intercept without eikonal-type corrections and the "bare" value, 
respectively. Both extremes are not excluded by experimental evidence The Pomeron slope has its 

usual value a' P = 0.25 GeV ~ 2 . The second function, the Pomeron structure function Fp, is proportional to 
the virtual photon-Pomeron cross section. It was emphasized in Ref. Q that other definitions of the "flux 
factor" are possible, e.g. differing from the one above by a constant. Therefore the normalization of the 
Pomeron structure function depends on the particular choice of the flux factor. For large values of Mx or 
small values of (3, the Mueller generalization of the optical theorem for inclusive cross sections can be used 
to represent the cross section given by Eq.(|l0|) in terms of the triple-Reggeon diagram shown in Fig. 5b. 
The Reggeon exchange in the upper part of the diagram is dominated by the Pomeron and / Regge poles. 
Since each term factorizes into a Reggeon propagator and a vertex, it is possible to obtain the structure 
function of the Pomeron from that of the deuteron Q] . 

We will below consider the diffractive production involving light quarks and heavy quarks separately, 
i.e. we write the Pomeron structure function as 

F P = F P + Ff . (15) 



For the contribution of light quarks to the structure function, we adopt the form for the deuteron structure 
function of Ref. |9j, but with the couplings of the exchanged Reggeon to the deuteron replaced by the 
couplings to the Pomeron, 

F P (P, Q 2 ,t) = e f d Cf((3, Q 2 )/3 1_Q/ (1 - f3) n ^~ 2 + e£c!(Q 2 )/r A(c?2) (l - /3)™W 2 )+ 2 , (16) 



with the ratios of the coupling constants 



k _ ' pp(t) [17i 



where r PP and g\ d are the couplings of the Pomeron (k = P) or the leading / Regge pole (k = /) to 
the Pomeron and to the deuteron, respectively. The values for the ratios e\ can be estimated from soft 
diffraction data. In our calculations we will use e d — e d — 0.07, which was shown in Ref. H to give a good 
description of the data on the diffractive structure function F% . The behaviour of the structure function 
for f3 — > 1 is determined by the exponent n(Q 2 ) which was obtained for the deuteron in Ref. |^] by using 
counting rules. For the Pomeron the exponent has to be adjusted accordingly since there is one parton 
spectator less. 

From experiment it follows that the triple Reggeon couplings are only weakly dependent on t and 
that for the Pomeron and the / Regge pole exchange this dependence is approximately equal. It can 



then be incorporated in that of the flux factor in Eq.(12), i.e. in g^ , which is then parametrized as 



9p P (t) — 9pp(®) exp(Ci) with C = 2.2 GeV~ 2 and g^,(0) 2 = 23 mb. With this parametrization, the structure 
function of the Pomeron (or its parton distributions) do not have the t-dependence any more. The values of 
the parameters mentioned above have been all taken from Ref. || . We again use the parametrization at an 
initial scale of Q 2 — 2 GeV 2 as starting point for QCD evolution to the desired Q 2 value. The contribution 
of the light quarks to the Pomeron structure function, Fp, is directly related to the corresponding singlet 
quark distribution of the Pomeron, 

F»(/J,Q 2 ) = ^ P (ftQ 2 ) , (18) 



with 

S P (/3,Q 2 )= E [/%M/3,Q 2 ) + /?g-M/?,Q 2 )] ■ (19) 
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This singlet distribution at Qq is shown in Fig. 6a. 

We now discuss the charmed quark contribution to the diffractive Pomeron structure function, which 
is related to the analogous charm quark distribution, 

F P s ((3,Q 2 ) = ^[(3c P ((3 : Q 2 )+(3c P ((3,Q 2 )] . (20) 

Using photon-gluon fusion as the leading mechanism, the diffractive charm production can be calculated in 
the same way as for charm production in DIS in Section 2. This mechanism for diffractive production is 
shown in Fig. 7a. Thus while charm production in DIS was a probe of the gluon distribution in the proton, 
the diffractive production of charm now probes the gluon distribution in the Pomeron, which is our main 
point of interest here. In addition to the "sea" contribution in Fig. 7a, there is also a "valence" contribution 
from the Pomeron structure function. A typical example of such a contribution can be seen in Fig. 7b, 
which has been studied in Ref. ]34j ]. The role of these additional mechanisms, which require further model 
assumptions, is not clear yet and had not been discussed in other papers on diffractive charm production. 
While the "sea" component is concentrated mostly at small values of /?, it can be shown that the "valence" 
contribution has its maximum at (3 ~ 0.5 in Fig. 7b. Its value is smaller by a factor of order 1/m 2 . In this 
region the "valence" contribution can therefore become comparable to the "sea" component. 

In the CKMT approach, the gluon distribution of the Pomeron can for low (3 be obtained from that 
of the proton (or deuteron), analogously to the Pomeron quark distribution discussed above. For a large 
parton momentum fraction, however, the analogy breaks down. The gluon distribution in the nucleon is 
determined at large x by the gluon radiation originating from the quarks; in the Pomeron, which is believed 
to be mainly composed of gluons, they are a priori present. In Ref. therefore the large (3 behaviour was 
modified according to 

/35p(/3,Q 2 ) = e^C 5 /3- A(Q2) (l-/3)" 9 , (21) 

where n g is a free parameter. The distributions we use are thus singular at j3 = due to the f3~ A dependence 
dictated by the Pomeron exchange. The observed Q 2 dependence of the data indicate that the gluon 
distribution should be rather hard and thus n g negative, i.e. between and —1 to yield a normalizable 
distribution. We use in the following two negative values of n g , and therefore our gluon distribution is 
also singular for (3=1. This can be seen in Fig. 6b, where we show our distributions for a scale of Q 2 = 2 
GeV 2 . For the actual calculation of diffractive charm production, we use QCD evolution to obtain the gluon 
distribution at the factorization scale of 4m 2 . 

Before showing our results for the diffractive production of light as well as charmed quarks, we first 
discuss the models of Refs. |6[j^] , where a similar description of the reaction involving the Pomeron structure 
function was used: model I of Gehrmann and Stirling (GS) and model 3 of Golec Biernat and Kwiecinski 
(GK) §. 

GK take the following singlet quark and gluon distributions, respectively, at Q 2 — 4 GeV 2 : 

H P (P) = 0.069 K /3 a44 (I - (3) 0M0 , (3g P {(3) = 1.16 K /3 5 . (22) 

The constant K relates the measured i-integrated diffractive structure function to the Pomeron structure 
function: 

F 2 D (x,Q 2 ,x P ) = J dt F 2 D (x,Q 2 ,xp,t) = K Q-^j F P ((3,Q 2 ). (23) 

The value for K depends on xp and is approximate 10 for the current experiments. The distributions used 
by GS at Q 2 , = 2 GeV 2 are: 

M/?) = 1-02 0(1-/3) , f3g P (P)= 4.92 /3(1 -/3) . (24) 



8 



In order to compare their quark and gluon distributions to ours, one has to take into account that different 
normalizations of the assumed Pomeron flux, f(xp,t), are used and they have to be converted to our 
conventions. 

In Fig. 6a and 6b we compare the different model distributions for quarks and gluons in the Pomeron, 
all normalized according to our definition of the flux factor. We show the distributions as originally 
parametrized on their respective initial scales, Q 2 ,. The quark singlet distribution, Sp, looks similar in 
the three models. Except for the behaviour of the CKMT parametrization at very low /?, the singlet quark 
distributions show roughly the same qualitative features for all models, because they were fitted to diffrac- 
tive production data for (3 > 0.065. Some differences with the CKMT model are due to the fact that the 
other models fit to all the diffractive data and thus also simulate doubly diffractive data with their fit, which 
account for about 30% of the cross section. 

In contrast to the quark distributions the gluon distributions of the three models are dramatically 
different: the predicted distributions have different shapes. The GK parametrization is the "hardest" of the 
three types of distributions, remaining essentially zero for low beta and having its main strength near (3=1. 
The distribution of GS is symmetric, vanishing at (3 = and (3=1, and thus is peaked in the middle. The 
gluons in the CKMT model are distributed with nearly constant density for intermediate (3 values; their 
distribution diverges at both ends. It becomes steeper at [3 = 1 the closer n g gets near —1. The very large 
values of gp((3) for intermediate values of (3 in the GS and GK models are difficult to reconcile with the 
cross sections for jet production in diffractive hadronic interactions |35| . The most direct test of all models 
is provided by the diffractive production of heavy quarks and we will therefore now look at diffractive charm 
production in DIS. 

For the calculations, the gluon distributions are used at the common factorization scale of 4mJ. They 
are shown in Fig. 6d and can be seen to have rather different qualitative features than in Fig. 6b. The QCD 
evolution from the initial scale Qq to the higher scale shifts the gluon distributions towards lower momentum 
fractions. This can be most clearly seen for the GS model. All gluon distributions are now singular at (3 = 0. 
For GS and GK this is due to the evolution, while CKMT already starts out singular at the origin due to 
the /3 _A behaviour, which is only somewhat modified by QCD evolution. The opposite effect can be seen 
at (3 = 1: the originally singular distribution in the CKMT model becomes zero. For completeness, we also 
show an example for the effect of evolution on the light quark distribution in Fig. 6c. For GS and GK, the 
magnitude of the distribution increases for moderate values of (3. Furthermore, again a slight shift towards 
small (3 occurs. For CKMT, the QCD evolution makes the density becomes more flat and the change in 
magnitude is such that the distribution is lower than the others at intermediate (3. 

One has to be careful when considering the energy-momentum sum rule, i.e. the integral over (3 of 
the sum of the gluon and quark distributions times (3. As was pointed out before, the choice of the Pomeron 
flux factor fixes the normalization of the Pomeron structure function and therefore the value of the energy- 
momentum sum rule. This sum rule needs not to be equal to 1 as was assumed in several publications. 
However, the value of this energy-momentum sum is preserved under QCD evolution, "virtual particle" . 
This explains the changes in magnitude discussed above: in response to an overall growth in the quark 
density the gluon density decreases. At the factorization scale, the total gluon content of the Pomeron 
predicted by CMKT is 60% and about 80% for the other models. 

The contribution of charm quarks to the Pomeron structure function, produced by the gluon-photon 
fusion mechanism is shown in Fig. 8 for the three gluonic distributions discussed above and for different 
values of Q 2 . The results are obtained by carrying out a convolution analogous to electroproduction. To 
obtain Fp c at a given momentum fraction (3 carried by the struck charmed quark, we have to integrate over 
gluon momentum fractions starting from a Q 2 dependent cut-off value that follows from the lower limit 
in Eq. ([!]). Both the CKMT and GS result increase continously as (3 decreases from this cut-off. The GK 
result has a maximum about halfway which originates from its relatively hard gluon distribution. As the 
virtuality of the photon increases, the cut-off moves to larger values, but the general shape of the curves 
doesn't change. For a Q 2 value as high as 500 GeV 2 , the factorization scale we have chosen may not be 
appropriate anymore and a Q 2 dependent value might be better. In general, the shown predictions are 
sufficiently different such that future experiments for charm production in hard diffractive scattering should 
be able to discriminate between them. 
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To see how significant the charm production contributes to the total diffractive cross section, we 
also show in Fig. 8 the full Fp. Adding these charm quark contributions results in a (small) Q 2 dependent 
violation of the sum rule discussed above. The fraction of the total structure function provided by the charm 
is in all cases most important for small j3 and increases with Q 2 . A measurement of the total diffractive 
cross section is of course a test of the different parton density input, contributions, 

The exponent n g , which determines the "hardness" of the CKMT gluon distribution, is a free pa- 
rameter. To see how sensitive the results are to its value, we compare in Fig. 9 results for Fp and Fp c for 
n g = —0.5 and —0.9 at the same Q 2 values as in Fig. 8. The charm part, Fp c , is most directly sensitive to 
the gluon distribution. The shape of the curves does not change much. However, as Fp c receives its con- 
tribution from gluons with high momentum fraction due to the presence of the cut-off in the convolution 
Eq. ([!]), the harder distribution, n g = —0.9, yields a significantly larger prediction. This is most pronounced 
at high P and similar for all values of Q 2 . For the total Fp, the more singular gluon distribution results 
in an increase, in particular at larger Q 2 . The exponent n g enters on the one hand through the charm 
contribution discussed above. On the other hand, it also influences the light quark contribution through 
the QCD evolution. The curves clearly show that the hardness of the gluon distribution is relevant for the 
total diffractive cross section. 

It is interesting to see the contribution of the diffractive production to the total DIS structure function, 
F2(x,Q 2 ). To compare to the results in Section 2, we convert this into the virtual photon cross section, 
<t(W, Q 2 ). For this purpose, we must integrate the diffractive contribution in Eq.(|l2]) over xp and t. For xp 
we integrate up to 0.1, which defines our lower limit for the rapidity gap. We use for the Pomeron intercept 
in the flux factor a value ap(0) = 1.13 as in Ref. ||. We show in Fig. 10 these diffractive contributions at 
the same Q 2 values as for the total in Fig. 4. They increase rapidly from their threshold and then reach an 
approximately powerlike VF-dependence. Comparison of Figs. 4 and 10 shows that the diffractive production 
is for large energies an order of magnitude smaller than the total. For large energies, the result for the harder 
gluon distribution is about a factor 1.5 higher. 



IV. BEAUTY PRODUCTION 



We now repeat some of the above calculations for the production of beauty quarks. Due to the larger 
mass, this of course involves a different factorization scale and thus probes the gluon distribution at a much 
higher scale; we take it as Am 2 . The cross section can be calculated using the same photon-gluon fusion 
mechanism and the approach discussed above. The only difference is the mass of the b-quark, which we 
take as m& = 4.7 GeV. It can be found using the mass of the charm quark from the following relation 

m& — m c — ffiB — mo , (25) 

which is valid in QCD up to small corrections ~ I/tuq p(i[ ]. The quantity ffiB = (tub + Stub* /4) in Eq.(p5|) 
denotes the center of gravity value for the lowest mesonic state (the same holds for mo). 

Predictions for the cross section for beauty production by real and virtual photons are shown in 
Fig. 11. The threshold is higher than for charm production, but the general shape of the cross section is 
similar. The magnitude is typically two orders of magnitude lower, which is due to the larger mass and the 
smaller charge. The contribution of diffractive beauty production to Fp(f],Q 2 ) is shown in Fig. 12. The 
(3 value for which the result becomes zero has decreased due to the larger quark mass. The shape of the 
curves, their dependence on Q 2 and the hardness of the gluon function are similar as for charm production, 
Fig. 9. The contribution to the structure function Fz{x, Q 2 ) is shown in Fig. 13. The threshold energy is now 
higher, but for the rest a similar behaviour to charm production can be seen: as in Fig. 10, the diffractive 
contribution is again an order of magnitude smaller than the total DIS beauty production and the sensitivity 
to the hardness of the gluon distribution looks the same. In general, these cross sections are much smaller 
than in the case of charm and large statistics is needed to observe them experimentally. 
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V. SUMMARY 



The main goal of this paper was to probe the gluon distribution in the Pomeron by means of diffractive 
heavy quark (charm, beauty) production. Given the present uncertainty in the data as well as the high 
sensitivity of perturbative QCD calculations to the heavy quark mass, we have here pursued a simple 
phcnomcnological approach to obtain a description of the main features of heavy quark production. First 
we showed that the data for total open charm photo- and electroproduction could be well described by 
using the gluon distribution of the proton predicted by the CKMT model, together with the photon-gluon 
fusion mechanism. From this comparison to the data we extracted a value of m c — 1.4 GeV for the 
charm quark mass parameter and confirmed the finding of Ref. |l5| that Am 2 is a good factorization scale. 
Having established this basis, we proceeded to diffractive charm production. We obtained results with gluon 
distributions of different models and compared their contribution to the total diffractive structure function 
of the Pomeron. For the CKMT model, we investigated how the results depend on the high momentum 
behaviour of the gluon distribution. A similar study was made for beauty production. 

In conclusion we found that different models for the gluonic content of the Pomeron lead to sizeable 
differences for diffractive charm production which experiments, such as the future high statistics experiments 
at HERA j37|, should be able to distinguish. For beauty production, similar features were found, but the 
cross sections are much smaller and thus more difficult to measure. 
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APPENDIX 

For the QCD evolution equation one needs the quark distributions at an initial scale Q\. Up to 
Q 2 = 5 GeV 2 the total valence and sea quark distribution in LO can directly be extracted from ^(x, Q 2 ) 
as parametrized by the CKMT model in Ref. ||: 

F™\ X ,Q 2 ) = ^{x,Q 2 ) + \d u {x,Q 2 ) = C f {x,Q 2 )x 1 -^{l-x) n{Q ^ , 

Fr a (x, Q 2 ) = lu s (x, Q 2 ) + ld s (x, Q 2 ) + ls s (x, Q 2 ) = C P (Q 2 ) X -^\l - x) n ^ +i , (26) 




This parametrization has been constructed in such a way that for x ~ 1 it is in accordance with the dual 
parton model at low Q 2 and with dimensional counting rules at very large Q 2 . The low x behaviour is readily 
explained in terms of Reggeon exchanges; the secondary Regge trajectory with intercept a/, corresponding 
to (/, j42)-exchanges, determines the small x distributions for the valence quarks, while the "effective" 
Pomeron exchange determines them for the sea quarks (and the gluon). The exponents in C/ and Cp are 
chosen such that the photolimit (Q 2 ,x — > 0) is finite. 
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The function B was written as the sum of the u and d quark contribution, B — B u + Bd, and is 
different for the proton and the deuteron: 



B U = \C U , B d = ~Q(1 - x) for the proton , 

B u = § [C u + C d {l -x)] , B d = i [C u + C d (l - x)} for the deuteron 



(28) 



where the constants C u and Cd were obtained from the condition that the integrals of the valence distri- 
butions over x should give the correct number of valence up and down quarks in the proton. Here the 
constants for the deuteron are related to those of the proton through isospin symmetry. For the sea quarks 
the simple assumption is used that u B = d s and s s — 0.5 u s , in reasonable accordance with the results from 
vN interactions. 

Also the gluon distribution for low x is determined by the Pomeron exchange and is thus proportional 
to the sea quark distributions. However, since sea quarks are produced mostly by gluons, the distribution 
of gluons will be harder than the one of the sea quarks. This leads to the gluon distribution in Eq.(Q): 

xg(x, Q 2 ) = G F \ (X > ° 2) = C g (x, Q 2 )x~^\l - x) n ^+ 3 , (29) 



where the proportionality factor G was obtained from the momentum sum rule at Q 2 — 2 GeV 2 . 

The parameter do in the effective Pomeron intercept was originally |9) put equal to 2 in order to 
have for the bare Pomeron a value of Ab are that is 3 times larger than A(0), in agreement with analyses 
of hadronic interactions (2^] . The seven remaining free parameters of the model were fitted in Ref. ||] to 
the data on total jp interaction cross section and to the proton structure function data of NMC. They are 
given by 

a = 0.2631 , 6 = 0.6452 , c = 3.5489 , di = 1.1170 , , . 

A = 0.1502 , A(0) = 0.07684 , a f = 0.4150 . ^ 

Here the parameters a, b, c and d\ are given in units of GeV 2 . With these values for the free parameters one 
finds explicitly for the coefficients: 

C u = 2.714 , C d = 1.618 , G = 12.63 . (31) 

However, using new HERA data in the region of Q 2 < 2 GeV 2 |38|,[59| it is possible to determine the 
parameters do and d\ with better accuracy: do = 2.2 , d\ — 0.6 GeV^. 
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FIGURE CAPTIONS 



Fig. 1 a. Photon-gluon fusion diagram. 

b. Charm production by the gluon-gluon fusion mechanism ("resolved" contribution). 

Fig. 2 Charm photoproduction cross section predicted by the CKMT model. Solid lines: direct contribution for 
different values of the charm mass; from top to bottom m c = 1.3, 1.4, 1.5, 1.6 GeV. Lower dashed line: resolved 
contribution with m c = 1.4 GeV obtained with the gluon distribution in the photon according to the CKMT model. 
Upper dashed line: the sum of the direct and resolved contribution (m c = 1.4GeV). 

Fig. 3 a. Charm contribution to the proton structure function F2(x,Q 2 ) in the CKMT approach as function of Q 2 
for different x; from top to bottom: x = 0.0001, 0.001, 0.01, 0.1. Solid lines are obtained by using the gluon fusion 
process below Q 2 = 50 GeV 2 and massless QCD evolution for Q 2 > 50 GeV 2 . Dashed lines: charm contribution 
by using the gluon fusion model for all Q 2 . Factorization scale is p, 2 — 4m 2 . 

b. Logarithmic derivative of obtained from Fig. 3a for x — 0.1 (dot-dashed), x = 0.01 (dotted), x — 0.001 
(dashed) and x = 0.0001 (solid line) 

c. Same figure as a., but with factorization scale fj, 2 = 4m 2 + Q 2 . 

d. Same figure as b., but with factorization scale /i 2 = 4m 2 + Q 2 . 

Fig. 4 Comparison of the charm photo-and electroproduction cross sections predicted in this paper with experiment 
for different Q 2 ; from top to bottom Q 2 = 0, 1.39, 2.47, 4.39, 7.81, 13.9, 24.7, 43.9, 78.1 GeV 2 . Solid lines: $ = 4m 2 . 
Dashed lines: ft 2 = 4m 2 + Q 2 . All curves are rescaled with powers of 10; the k-th curve from the top is rescaled 
with a factor 10~ fe . The data are from |||^ofl3|] . 

Fig. 5 a. Diffractive dissociaton of the photon in photon-proton scattering due to the Pomeron exchange, 
b. Triple-regge diagram for hard diffractive scattering in the small (3 region. 

Fig. 6 a. Different parametrizations for the light quark distribution of the Pomeron, Ep(/3), at initial scales Q 2 ,. 
Solid line: this paper with n g — —0.5. Dotted line: this paper with n g = —0.9. Dot-dashed line: GS Q. Dashed 
line: GK §. 

b. Different gluonic distributions in the Pomeron at initial scales Q 2 ,. Solid line: this paper with n g = —0.5. Dotted 
line: this paper with n g — —0.9. Dot-dashed line: GS Q. Dashed line: GK [Q. 

c. Predictions for the light quark singlet distributions at the factorization scale by different models. Solid line: this 
paper with n g — —0.5. Dotted line: this paper with n g — —0.9. Dot-dashed line: GS Q. Dashed line: GK |J. 

d. Predictions for gluon distributions at the factorization scale by different models. Solid line: this paper with 
n g = —0.5. Dotted line: this paper with n g — —0.9. Dot-dashed line: GS Q. Dashed line: GK jS. 

Fig. 7 a. Pomeron "sea" charm contribution. 

b. Typical diagram that contributes to the "valence" charm component of the Pomeron. 

Fig. 8 Pomeron structure function (upper curves) and its charm contribution (lower curves) predicted by different 
models as function of (3 at different Q 2 . Solid line: this paper with n g = —0.5. Dot-dashed line: GS H. Dashed 
line: GK §. 

a. Q 2 = 10 GeV 2 . 

b. Q 2 = 100 GeV 2 . 

c. Q 2 = 500 GeV 2 . 

Fig. 9 CKMT predictions with different n g for the Pomeron structure function (upper curves) and its charm 
contribution (lower curves). Solid lines: n g = —0.5. Dotted lines with n g = —0.9. For every set the Q 2 -values are 
from top to bottom: Q 2 = 500, 100, 10 GeV 2 . 

Fig. 10 Diffractive charm contribution to the total charm production cross section for different values of Q 2 ; from 
top to bottom Q 2 = 0,1.39,2.47,4.39,7.81,13.9,24.7,43.9,78.1 GeV 2 . Solid lines: n g = -0.5. Dotted lines: 
n g = -0.9. 

Fig. 11 Cross sections for photo-and electroproduction of beauty predicted by this paper for different values of Q 2 ; 
from top to bottom Q 2 = 0, 1.39, 2.47, 4.39, 7.81, 13.9, 24.7, 43.9, 78.1 GeV 2 . 

Fig. 12 Beauty contribution to the Pomeron structure function predicted by this paper at different values of Q 2 : 
from top to bottom Q 2 = 500, 100, 10 GeV 2 . Solid lines: n g = -0.5. Dotted lines: n g = -0.9. 

Fig. 13 Diffractive beauty contribution to the total beauty production cross section predicted by this paper for 
different values of Q 2 ; from top to bottom Q 2 = 0,1.39,2.47,4.39,7.81,13.9,24.7,43.9,78.1 GeV 2 . Solid lines: 
n g = —0.5. Dotted lines: n g — —0.9. 
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